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' We discuss the onset of symmetry breaking from the false vacuum in generic scenarios in which 

| the mass squared of the symmetry breaking (Higgs) field depends linearly with time, as it occurs, 

via the evolution of the inflaton, in models of hybrid inflation. We show that the Higgs fluctuations 

■ evolve from quantum to classical during the initial stages. This justifies the subsequent use of real- 
time lattice simulations to describe the fully non-perturbative and non-linear process of symmetry 
breaking. The early distribution of the Higgs field is that of a smooth classical gaussian random field, 
and consists of lumps whose shape and distribution is well understood analytically. The lumps grow 

^ with time and develop into "bubbles" which eventually collide among themselves, thus populating 

the high momentum modes, in their way towards thermalization at the true vacuum. With the help 

■ of some approximations we are able to provide a quasi-analytic understanding of this process. 
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The problem of symmetry breaking in quantum field theory (QFT) has been with us for several decades. In 
the context of cosmology, it has usually been associated with thermal phase transitions and the production of 
topological defects |2|. Understanding the way in which the order parameter associated with the breaking of the 
CU symmetry evolves from a symmetric state (the false vacuum) to a broken state (the true vacuum) is non-trivial 0, 0, 
LJilIl] ■ Only recently has this problem been addressed in the context of symmetry breaking at zero temperature, at the 
| end of a period of hybrid inflation B, |p , in the so-called process of tachyonic preheating, i.e. spinodal instability, in the 
context of preheating after inflation [lOj. There, classical evolution equations have been solved with real-time lattice 
simulations, developed for studying the problem of preheating [ill Il2t IT^|. which include all the non-perturbative 
• i-H , and non-linear character of the phase transition. It was found that symmetry breaking occurs typically in just one 
oscillation around the true vacuum 0], most of the false vacuum energy going into gradient modes, rather than 
kinetic energy. 

However, the problem of the transition from a false quantum vacuum state at zero temperature (as occurs at the end 
of a period of inflation) to the true quantum vacuum state full of radiation at a certain temperature has not been fully 
addressed yet. Most of the previous approaches refer to the decay from a false vacuum state at finite temperature, 
relying on the Hartree or large N approximation 0, 0, 0, 0- Only recently, the zero temperature problem was 
addressed within the classical approximation but mainly for the case of an instantaneous quench [T^ | , which can lead 
to cosmologically interesting particle production [l5| . Tachyonic preheating was recently studied beyond the quench 
approximation in Ref . fl6| . 

We will argue that symmetry breaking proceeds through a state in which the relevant degrees of freedom are semi- 
classical infrared modes, which can be described in a non-perturbative and non-linear way with a classical effective 
field theory, whose classical equations of motion can be solved numerically in the lattice and thus allow us to study 
the fully non-perturbative out of equilibrium process of symmetry breaking. 

The quantum to classical transitions of field Fourier modes have been addressed before in the context of inflation 
[r?l llH , EH l2fj| , where it is mandatory to understand the transition from quantum fluctuations of the inflaton field 
during inflation to the classical metric fluctuations on super-horizon scales, since they are believed to be responsible for 
the observed temperature anisotropies in the cosmic microwave background, as well as the scalar density perturbations 
giving rise to galaxies and large scale structure formation. The use of the classical approximation to study the process 
of preheating after inflation has been proposed in and, in a context similar to ours, it has been recently used in 
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In this paper we will use such a well developed formalism to study the first instances of a generic symmetry breaking 
process, i.e. the conversion of quantum modes of the symmetry breaking field (the QFT order parameter, generically 
called the Higgs field) into a classical gaussian random field whose subsequent non-linear evolution equations can be 
solved with lattice simulations. We give here a self-contained presentation of the conditions under which classical 
behaviour holds, and apply it to the analysis of the false vacuum decay after inflation, with specific initial conditions. 

Our final aim is to study the process of electroweak symmetry breaking, and the possibility of realizing baryogenesis 
at the electroweak scale, via the non-equilibrium process of preheating after inflation [2lLl22ll2^ |. Therefore, this paper 
is intended as the first one in a series, in which we will progressively incorporate more complexity, i.e. gauge fields, 
Chern-Simons, CP-violation, etc., into the picture. Of course, our results are readily generalisable to any other phase 
transition that may have occurred in the early universe at the end of a period of hybrid inflation, e.g. at GUT scales. 

The paper is organised as follows. In Section II we describe the initial conditions for spontaneous symmetry 
breaking coming from a hybrid model of inflation. The inflaton acts here like a background field whose coupling gives 
a time-dependent mass to the Higgs. In Section III we study the quantum evolution of the Higgs field in the linear 
approximation from the bifurcation point. The Fourier modes decouple in this approximation and can be studied 
as a quantum mechanical ensemble of harmonic oscillators, both in the Heisenberg and the Schrodinger picture. We 
then study, with the use of the Wigner function, the quantum to classical transition of the Higgs modes. We show 
that each quantum mode can be described exactly like a classical gaussian random field, and give a prescription for 
computing the Weyl-ordered quantum expectation values of operators in terms of classical averages over a gaussian 
random field with the Wigner function as probability distribution. We then define and characterise when a mode can 
be called "quasi-classical". 

In Section IV we give the exact solutions to the field evolution equations of the Higgs in the linear approximation in 
terms of Airy functions, and show that soon after the bifurcation the infrared modes become quasi-classical according 
to the definition of the previous section. This analysis follows closely, although in greater detail, what has been 
studied previously in the literature |l6l l33l l34j . In Section V we analyse the inclusion of the non- linear terms in the 
quantum evolution within perturbation theory. We give a prescription for treating the ultraviolet divergences and to 
renormalise the parameters of the theory. This leads to a regular probability distribution to be used for a classical 
field description, which matches the renormalised quantum expectation values of the Weyl ordered products. The 
matching is done at a time in which the infrared modes have grown sufficiently to be well described as classical modes. 
This occurs well before non-linearities are important, and therefore our gaussian approximation is valid, in a similar 
spirit as that of Ref. \J\. The quantum ultraviolet modes, on the contrary, can be thought as integrated out, and used 
to renormalise the parameters of the classical theory. We might then interpret our classical field distribution as an 
effective theory for the long wave-length modes. We also estimate the time at which symmetry breaking sets in. 

In Section VI we describe the methodology to be used to take care of the full non-linear evolution of the system. 
The initial space-time structure of the classical Higgs field is analysed in VI. A. Being a gaussian random field, it can be 
described in a similar way to the matter density field whose fluctuations give rise to galaxies and large scale structure 
via gravitational collapse (24[. The Higgs field is found to possess, at symmetry breaking, an inhomogeneous spatial 
distribution made of lumps, whose shape and initial evolution can be well understood analytically. The space-time 
inhomogeneous character of symmetry breaking in the Higgs-inflaton system has also been reported in Ref. |16| for a 
one component Higgs model. For the complete non-linear dynamics of the full Higgs-inflaton system we make use of 
lattice real-time evolution methods. The details of our procedure, its connection and difference with lattice methods 
used by other authors, and a detailed check of the validity of the approximations used, are described in Section VLB. 

Finally, in Section VII we present the results of this stage of the evolution of the system. The lumps mentioned 
in the previous sections grow and, once its center reaches the Higgs vacuum expectation value, invaginate and create 
an approximately spherically symmetric "bubble" which expands at a very high speed. Meanwhile, the center of the 
lump-bubble continues to oscillate with decreasing amplitude, leading to secondary bubbles. All these phenomena 
can be well understood with the help of some approximations which reduce the full non-linear equations to a one 
or two dimensional partial differential equation of a single scalar field. This simplified picture matches qualitatively 
and (to a high degree) quantitatively the results of the lattice simulations. Eventually, bubbles centred at different 
points collide and transfer most of their potential and kinetic energy to gradient energy, thus populating the higher 
momentum modes. This process leads to complete symmetry breaking and (classical) thermalization. The whole 
history of the system is illustrated by following the evolution of 2D sections of a particular configuration. We also 
show histograms for the field values of both Higgs and inflaton, which start in the false vacuum and are seen to end 
up peaked around the true vacuum. In Section VIII we draw our conclusions and describe the future directions in 
which this work can be extended, first by including the production of SU(2) gauge fields and afterwards by studying 
the rate of sphaleron transitions that may give rise to a non-negligible amount of baryons. 

We have added three appendices. In Appendix A we describe the formalism of squeezed states following [T§. fli^ . 
which can be applied to the initial stages of the Higgs evolution in the linear regime, and gives rise to the semi-classical 
nature of the long wavelength modes. In Appendix B we compute the Wigner function for the evolved gaussian initial 
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vacuum state |l8l Il9l l20j . and show explicitly the squeezing of the infrared modes. We give a definite condition for 
characterising the moment in which a mode can be treated as quasi-classical. In Appendix C we give the details of 
the perturbative calculations of the non-linear evolution of our system, both at a classical and quantum-mechanical 
level. 



II. THE HIGGS FIELD AT THE END OF HYBRID INFLATION 



The precise model of hybrid inflation will not be important for our purposes here. However, for concreteness we will 
implement it in the context of a super-symmetric extension of the standard model in which the radiative corrections 
are responsible for the running of the inflaton field during the few (~ 5 — 10) e-folds necessary to cool the universe so 
that the electroweak symmetry breaking (EWSB) occurs at zero temperature. The fluctuations responsible for CMB 
temperature anisotropies and large scale structure come from a previous stage of inflation, completely independent 
of this. Moreover, since EWSB occurs at low energies, we can, and will in what follows, safely ignore the rate of 
expansion, H ~ 1CP 5 eV, during symmetry breaking and treat the fields as if they were in Minkowski space. In 
particular, one can consider the super-symmetric hybrid model of Dvali, Shaft and Schaefer |25j . where the super- 
potential fixes a relation between the couplings, g 2 = 2A. As we will see, this choice simplifies some stages of the 
dynamics of symmetry breaking after inflation |26j ] , but is not crucial. A study of the process of tachyonic preheating 
after a variety of more general super-symmetric models of inflation will be given in Ref. |27| . 

The hybrid model we are considering is a simple generalisation of the Standard Model symmetry breaking sector, 
which consists of the Higgs field, $ = ^(<j>o U + i<j> a T a ), with r a the Pauli matrices, and an inflaton x, a singlet under 
SU(2). The inflaton couples only to the Higgs, with coupling constant g. The scalar potential has the usual Higgs 
term plus a coupling to a massive inflaton, 

c = (D^y + l -(d^ x ) 2 - a(V $ - y) 2 - sV* f * - ~mV , (i) 

where v = 246 GeV is the expectation value of the Higgs in the true vacuum, fj, is the mass of the inflaton in the 
false vacuum and m = \f~\v. We are assuming implicitly that whenever there is a contraction O^O, we should 
take the trace over the SU(2) matrices, i.e. $t$ = Tr$t$ = 1(^2 _|_ a a ) = \cf)\ 2 /2. The Higgs mass in the true 

vacuum is determined by its self-coupling: mjj = \/2Xv, while the mass of the inflaton in the true vacuum is given 
by mi = > //. 

In this paper we will simplify the analysis of the dynamics by omitting the SU(2) gauge field and working with a 
generic Higgs field with iV c real components. We anticipate that the most important conclusions of this paper are 
not affected by the introduction of the gauge field, and leave for a forthcoming publication the symmetry breaking 
dynamics in the presence of gauge fields. For ease of notation we will drop the internal indices of the Higgs field 
whenever all components behave in the same way. The numerical simulations that will be presented correspond to a 
N c — 4 component Higgs field. 

During hybrid inflation |8( the Higgs field has a large and positive effective mass squared due to its coupling to the 
inflaton field, which slow-rolls down its potential valley. The potential for the coupled fields is the following 

vfax) = \{W - + V + ■ (2) 

where the parameters in the potential depend on the number of Higgs components N c as: A = Xo/N c , g 2 = g^/N c 
and v 2 = N c v 2 , with this m 2 = Xv 2 = \§v 2 is independent of the number of components. 

It is the effective false vacuum energy Vo = Xv 4 /4 = m 2 v 2 /4 which drives the period of hybrid inflation. Inflation 
ends when the inflaton homogeneous mode, x = (x); slow-rolls below the bifurcation point x — Xc = Tn/g, at which 
the Higgs is massless, 

Below the critical point, the Higgs has a negative mass squared and long wave modes will grow exponentially, driving 
the process of symmetry breaking [l4| . The process by which the mass squared of the Higgs goes from large and positive 
to large and negative is not instantaneous, but depends strongly on the velocity of the inflaton at the bifurcation 
point, 
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Typically the speed of the inflaton is such that the process takes place in less than one Hubble time, a condition known 
as the "waterfall" condition 0, [t| , which ensures the absence of a second period of inflation after the bifurcation 
point p8j . The actual value of V depends very much on the model and the scale of inflation, and we will treat it 
here as an arbitrary model parameter. In this case, the effective mass of the Higgs across the bifurcation point can 
be written as a time-dependent mass 

m%(t) = -2Vm 3 (t - tc) + 0(V 2 (t - t c ) 2 ) . (5) 

Note that a similar situation arises in the case of simple extensions of the standard model Higgs, in which radiative 
corrections (dominated by the large top quark Yukawa coupling) induce the running of the Higgs mass square from 
positive to negative thus providing a mechanism for Electroweak symmetry breaking. The role of the running scale is 
played here by the inflaton homogeneous mode. Alternatively, one can envisage a secondary period of hybrid thermal 
inflation [2|J,[3(j just above the electroweak scale, which lasted only a few e-folds and super-cooled the false vacuum, 
leaving only the fast rolling inflaton coupled to the Higgs. This short second period of inflation would not affect the 
CMB anisotropies, but would provide a natural initial condition for the growth of quantum fluctuation of the Higgs 
field, as they evolve across the bifurcation point, toward symmetry breaking. 

Let us consider the effective action for the Higgs field <I>(x, t) ignoring the self-coupling A-term (we omit the internal 
indices of the Higgs field), 



5 / d 3 xdt - 



(4>) 2 -m 2 -ml(t)<p 2 , (6) 



where we have included the time-dependent mass |JS}, to linear order, which is the only effect that the presence of the 
homogeneous mode x(i) of the inflaton field induces in the evolution of Higgs quantum modes. 
We now define a new scale M = (2V) 1 / 3 m, and thus redefine our coordinates as 

r = M{t-t c ) — ► = M^, (7) 

X = Afx — ► K = — . (8) 

M v ; 

where primes denote derivatives w.r.t. r, and k is the wavenumber associated with the Higgs Fourier modes, 

d 3 x 



f d 3 x 

$(k,r) = J /2 $(x,t) exp(-ix-k). 



From now on, we will use x and k as the normalised position and momentum coordinates, i.e. we will work in units 
of M = 1. We will also denote the normalised Higgs quantum fluctuations by y = <j)/M, for which the effective action 
is 

5 = |d 3 xdri[( 2/ ') 2 -(Vy) 2 + r 2/ 2 ] . (9) 
dC , 

We can define the conjugate momentum as p = — — = y , and thus the corresponding Hamiltonian becomes 

dy' 

H = j >xi [ J3 2 + (V 2/ ) 2 -ry 2 ] . (10) 
In momentum space, the Hamiltonian becomes 

H = J d 3 ki [p(k,T)pt(k,r) + (fc 2 -r)y(k,r)yt( k , T )] . (11) 

The Euler-Lagrange equations for this field can be written in terms of the momentum eigenmodes as a series of 
uncoupled oscillator equations: 

y"(k,r) + (fc 2 -r) 2/ (k,r) =0. (12) 



III. QUANTUM EVOLUTION IN THE GAUSSIAN APPROXIMATION 

In this section we will start the description of the quantum evolution of the system assuming that we can neglect 
the non-linear terms which are proportional to A. Our goal is to determine the precise conditions under which the 
system evolves into a classical one. Our presentation will be general and applicable to any time dependent harmonic 
oscillator system with time-dependent spring constant w 2 (fc, r), only in the next section we will apply this formalism 
to our particular problem (w 2 (fc, r) = k 2 — r). Our results overlap and coincide with Refs. |l7llla.ll9|. 
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A. The Heisenberg picture 



In the Heisenberg picture the quantum system is described by means of the position y(k, r) and momentum operators 
p(k, r) corresponding to each oscillator. The canonical equal-time commutation relations for the fields (h = 1 here 
and throughout) in position and momentum space are 



y(x,r), p(x',r) = iS 3 {x~x'), 



y(k,r), p(k',r) = iS 3 (k + k') 



(13) 



Furthermore, hermiticity of the operators in position space imply the relations y^(k, r) = y(— k, r) and p^(k, t) = 
p(-k,r) 

We will assume that at r = r = 0, i.e. at the bifurcation point t = t c , the state of the system is given by the 
ground state of the Hamiltonian with oscillator frequency w(fc,0) — k. It is then useful to express the position and 
momentum operators in terms of creation-annihilation operators at that time: 



2/(k,T ) = -j= (a(k,T ) +a t (-k,r )) , 
p(k,T ) = -i (a(k,T ) - a f (-k,T )) 



(14) 



The quantum operators satisfy the classical equations of motion, which we will write down as a system of coupled 
first-order equations 



dr V ' T ~ dr \y(k,T) 
whose solution can be expressed as: 

v(k,r) = M(k, T)v(k,T ) = 



-c 2 (fc,r)\ /p(k,r)' 

1 J L(k,r) 



(15) 



fffci(r) V2kg k2 (T) 



k fk2[r) V2kf kl (T) 



v(k,r ) . 



(16) 



where f k i = Re/fc and = Im/fc, with /fc(r) a complex solution of the equation of motion, with initial conditions, 

fl + r)) / fc = , / fe (r ) = , (17) 



and 



9k = 9k\ + igk.2 =if k - 



9k(T ) = \j - . 



(18) 



Note that since the motion is Hamiltonian (i.e. canonical), the determinant of M(fc, r) is det M(fc, r) = 1 ,Vr, a 
condition that is equivalent to the Wronskian of Eq. (|17fl being one at all times, 



* (f'k n - f'k'fk) = 9k ft +g* k f k = 2Re (g k f* k ) = 1 



(19) 



The previous formulae allow us to compute the expectation value of products of fields at any time r in terms of the 
expectation values of fields at time tq. Substituting l|14fl into the expression for the fields at time r we obtain 



y(k, t) = f k (r) a(k, r ) + f* k (r) at(-k, r ) , 
P( k . r ) = (3fc( T ) a ( k ,' r o) - .9fe(r)a t (-k,ro)^ 



(20) 



The quantum information of the system is encoded in the expectation values of products of fields. For a gaussian 
field the only quantities needed to describe the system are the two-point expectation values, 



(0, T k(k, r) v b (k', r')|0, r ) = E ah (fc, r, /) S 3 (k + k') , 



(21) 
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where |0, To) is the initial vacuum state satisfying a(k, To) |0, To) = , Vk. The value of this matrix at any pair of times 
can be expressed in terms of the matrix M and the corresponding expectation values at time t as follows: 

E(k,T,T')=M(k, T)H(k,T ,T Q )M T (k, T'). (22) 

The quantum initial condition on the state of the system at time To amounts to: 

s(fc,T ,T ) =h 2 y ^ 

Note that this matrix is hermitian, but neither real nor symmetric, and its determinant vanishes. The imaginary part 
results from the equal time commutation relations and does not depend on the particular state of the system. The 
real symmetric part alone characterises completely the state. 

Let us conclude this section by giving the expression of the equal time expectation values at any other time: 

/ | 5fe (T)| 2 F k (r)-i\ ( |. 9fc (T)| 2 _i^(T)|/ fe (T)| 2 \ 
£fc,T,T)= = (24) 

\F k {r) + i |/ fe (T)| 2 J \in k {T)\f k {T)\ 2 |A(t)| 2 J 

with 

g* k (r) l-2zF fc (T) 

nk(T) = m = ^mr> (25) 

F k (r) = Im(ttg k ). (26) 

As a consequence of the unit determinant of M(fc, t) one concludes that the determinant of the symmetric (real) part 
of E(fc,T, t) is time-independent and equal to 1/4. Note that using (|25(l we can rewrite the conjugate momentum as 



p(k,T)=p(k,T) + r ^^ 2 ,(k,T) (27) 



with 

p(k,r) = - \f k L\\2 (//c(r)a(k,T ) - f k (t) a f (-k, t )) , 
a relation that will prove useful in the next section. 



B. The Schrodinger picture and the classical limit 



Let us go now from the Heisenberg to the Schro ding er representation, and compute the initial state vacuum 
eigenfunction ^^(t = To). We will follow here Refs. [171 Ha. fl9|. In what follows we will denote operators in the 



Schrodinger representation by 
condition 



y(k, To) and = p(k, tq). The initial vacuum state |0, tq) is defined through the 



Vl + T 



Vk, a(k,T O )|0,T O ) = 

1 d 



key? 



^o{yly k \r Q ) =0 



2 m + l Wk 



Pk 



\0,T ) =0, 



(28) 



where we have used the position representation, j/k 



- , and Nq gives the corresponding normalisation. 



We will now study the time evolution of this initial wave function using the unitary evolution operator U — U(t, To), 
satisfying U' = —iJiU. The state evolves in the Schrodinger picture as |0, t) = U\0,tq). We can make use of the 
result of the previous section to determine this state. By inverting (|20[1 we find 



a(k, tq) = g* k (t) y (k, t) + i f k (r) p(k, t) , 



(29) 
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which acting on the initial state becomes, Vk , Vr, 

fk(r) 



U 



9t(r) 



p(k,r) 



to) 



10, r) 



0* 



K\fk(r)\ 



-nUr)\y° k \ 2 



(30) 



with r2fc(r) given by Eq. 1)25(1 . We see that the unitary evolution preserves the gaussian form of the wave functional. 
The wave function i|3U|) is called a 2-mode squeezed state. The normalised probability distribution, for each mode k, 



'o(«V 



1 



H/ fc (r)| ; 



exp 



( \vl\ 2 \ 



\fk(rW 



(31) 



is a gaussian distribution, with dispersion given by |/fc| 2 . This agrees with the result obtained in the previous section 
in the Heisenberg picture. The phase a(fc, t) cannot be determined by this method, but as we have seen it has no 
effect on the probability distribution nor on the Wigner function, see below. However, from the Schrodinger equation, 
idr^oir) = Wo(t), one can deduce that a'(k, r) = (2|/ fc (r)| 2 )- 1 . 
We can also compute the occupation number, n k , 



n k (r) = (0,r|at(k,r o )a(k,ro)|0,r) = ^\g k \ 2 + ^\f k \ 2 - \ 



(32) 



a quantity that is always positive definite. 

We now address the problem of approximating the quantum evolution just described by a classical evolution. For 
that purpose the vacuum expectation values of products of position and momentum operators should be recovered as 
ensemble averages of random fields. It is clear that for the non-interacting theory (A = 0) that we are considering, 
such a classical random field should be gaussian, with all the information encoded in the real expectation values 
of products of two fields. Only the symmetrical part of S(/c,t, t) is real, see (I24|l . and thus a natural candidate 
to be approximated by the classical gaussian random field. Notice that this corresponds to matching Weyl-ordered 
(symmetrised in y k , pi) quantum expectation values of operators through, in the Schrodinger picture, 



<0,r|G(y k ,p k )|0,r) w = (G(y k ,p k )), 



(33) 



where (0, r|G(y k ,p k )|0, t)w denotes the quantum average of the Weyl-ordered operator in the state given by the wave 
function l|3UI) . and (G(y k ,p k )) gs denotes the classical gaussian average. The latter is obtained as an average over a 

( FkW) \ 

gaussian ensemble with y k and p k — p k — y k independent gaussian variables with probability distribution 

V \fk{T)¥_J 

given by the Wigner function in phase space, see [18|, [ijj , equation lj2T|) and the Appendix B 



W ok (y,p) 



exp 



\fk 



A very trivial illustration of this equality is given by the following symmetrised vacuum expectation value 



2 (0, t\ £kp k / + p k yu |0, r) = F k (r)S 3 (k + k') 



while 



!yk (pk> 



F k (r) v. 
\h? Vk ' ))sS 



F k (r) 

IA-1 2 



y k 2/k') gs =F fc (r)<S 3 (k + k') 



(34) 



(35) 



(36) 



Even though only the symmetrised expectation values are described by the gaussian ensemble average, for the 
gaussian ground state of (|30|l quantum expectation values with arbitrary ordering of operators can also be computed. 
To be specific, the expectation value of any operator G(y,p), with any given ordering of y and p, can be rewritten 
as a lineal combination of Weyl-ordered operators with coefficients proportional to the commutator which is a time 
independent c-number; schematically 



(0, T\G(y,p)\Q, t) - (G (y, P )) gs + W (G n (y,p)) 



gs ■ 



(37) 
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where we have introduced fi as an expansion parameter to make explicit the connection with the semi-classical 
approximation. For instance, for the example in l|35[l we would obtain 

(0,T|y k p k ,|0,r) = (V fe (r)-£)<5 3 (k + k') (38) 

In this spirit, a quasi-classical state can be defined as a state for which the leading term in equation (|37|) dominates, 
and quantum averages can be approximated by 

(0,r\G(y,p)\0,r) « (G (y,p) ) gs . (39) 

This generically happens when (0, r\py |0,r)v7 ^> | (0, t| [p,y] |0,t)|, i.e. when the so-called WKB phase, Fk(r) in 
(|26|l . verifies |Ffc(r)| ^S> 1. For such a quasi-classical state the ambiguity in the ordering of operators is quantitatively 
negligible and classicality in the sense of (|39|) holds. As an illustration let us compute the following expectation value 
(for ease of notation we have omitted the k dependence of the operators and the delta functions) 

l(0,r\p 2 f + y 2 p 2 \0,r)^3F 2 (r)-^. (40) 

In the classical approximation we would obtain 

(p 2 y 2 ) gs = (P 2 ) es (y 2 ) g s + 2{py)l s = |/ fe (r)| 2 | 5fe (r)| 2 + 2F 2 (r) = 3F 2 (r) + ± (41) 

which reproduces the Weyl ordered part of the quantum result and is a very good approximation as long as | (t) j 3> 1. 

This approach works as long as the theory is non-interacting and the gaussianity of the quantum state is preserved 
by the evolution. In the problem at hand, we can assume this to be the case in the first stages of evolution before 
non linearities have set in, but not when the self-coupling term starts to be relevant just before symmetry breaking. 
However, as long as l-F/^r)! ^ 1, the state can still be approximately described, through the non-linear stages, via 
a classical random field. We will describe below in detail how this classicality follows for our specific problem (see 
also the Appendices A and B for the precise formulation of the squeezed states and the Wigner function formalism, 
following Refs. [H E HJ ) . 



IV. EXACT SOLUTIONS OF THE FIELD EVOLUTION EQUATIONS 



Let us apply now the above formalism to the case of the quantum fluctuation modes of the Higgs at symmetry 
breaking. This case was studied previously in [la, H2> liiil lii-l • For oj 2 (A,t) = fc 2 — t the linear equation (fT7| for the 
quantum modes of the Higgs field becomes 



fk+(k 2 -r)f k =0, with / fe (r = 0) 

Its solution can be given in terms of Airy functions |3l| : 



1 



2k 



A(t) = d(fc)Bi(r - fc 2 ) + C a (fc) Ai(r - fc 2 ) , 
Sfc(r) = iC 1 (fc)Bi'(r-fc 2 ) + iC 2 (fc)Ai'(r-fc 2 ), 



Cx{k) 
C 2 (k) 



'2k 



Ai'(-fc 2 ) + ifcAi(-fc 2 ) 



Bi'(-fc 2 ) +ifcBi(-fc 2 ) 



(42) 

(43) 
(44) 
(45) 
(46) 



which satisfy the Wronskian condition, g k + gl fk = 1, where we have used the corresponding Wronskian for the 
Airy functions, 

7T [Ai(z)Bi'(z) - Bi(z)Ai'(z)] = 1 . (47) 
We can then compute the occupation number, n k , (|32|l and the imaginary part of the WKB phase, i 7 /., see Ij26(l . 
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A. Long wavelength quasi-classical modes 

Still remains to be computed the time at which fluctuations become classical in the sense of Eq. (|39ll . As we will 
see, the field fluctuation modes will become quasi-classical as their wavelength becomes larger than the only physical 
scale in the problem, the time-dependent Higgs mass, i.e. A = 27r/fc 3> litjyfr. In order to show this, let us take the 
limit i; 2 «T for the long wavelength modes in the exact solutions l|43|) and (|44|l . 

fk(r) = Ci(fc)Bi(r) + C 2 (fc)Ai(r) ~ C t {k) w(r), (48) 
fffe (r) = id{k)B\!{r) + iC 2 {k)M{T) ~ iC\{k) w'(t) , (49) 

where the function w(t) — Bi(r) + \/3 Ai(r) is the one appearing in the Appendix A, and we have used the fact that 
C%{k) = \fiCi{k) in the limit k — > 0. Using the large z ^ 1 approximation of the Airy functions, see Ref. 31], 



1 



Bi( z )~ -== z'^e+i* 3 ' 2 , Ai^^ir 1 /^-!^, (50) 



'7T 



2V^F ' 



Bi'(z) - z 1 / 2 Bi(z) , Ai'(z) ~ -z 1 / 2 Ai(z) , (51) 

we conclude that the first terms in both //- and g k correspond to the growing modes, while the second terms are the 
decaying modes, and can be ignored soon after the bifurcation point. 

We are now prepared to answer the question of classicality of the modes. The wave function phase shift is given by 

F k = ImtffcSk) * \Cx{k)\ 2 w'{r)w{r) ~ | /fc |2 ~ | Cl (fc)| 2 I , (52) 

«;(TJ 7r 

which grows faster than exponentially at large time. On the other hand, the occupation number (|32|) is 

2k(n k + 1) = | 5fe | 2 + fc 2 |/ fc | 2 ~ (^) 2 |M 2 ~ ^ |d(fc)| 2 i e^ 3/2 . (53) 



Therefore, we have 



\F k (r)\d^±»l — > £ t (r),h t (r)»Kr)|, (54) 
(w'/w) 

that is, |Ffc| ^> 1 whenever the energy of the mode E k is much greater than the Higgs mass, computed as the 
instantaneous curvature of the Higgs potential. 

Notice in particular that under the condition \F k \ 3> 1, the momentum and field eigenmodes are related by 

9k{r) =n|(r)/ fe (r) — ► i r^rp /fc( T ) ■ 

In terms of the gaussian random fields, the momentum distribution of the Wigner function becomes a delta function, 

/ Fk (r) \ 
Hp~ mo- y i see Appendix B. 

v l/fc| 2 y 

We show in Fig. [T] the exponential growth of the phase Fk(r) as a function of momenta k, for different times. These 
plots were obtained using the exact Airy function solutions. Note that in the limit of large wavelengths fc 2 <C r, it is 
indeed verified that |Ffc(r)| ~ 2kn k (T)(w/w') 3> 1, as stated above. 

We can now compute the time for which a single mode k becomes quasi-classical, in the sense (|54f) . We have 
confirmed that after r ~ 2 modes with < fc < ^/r , which, as we will see later, is the range of interest, have become 
quasi-classical. We have drawn the line separating classical from quantum modes in Fig. |2 as a function of the mode 
fc. The high energy part of the spectrum always remains in the quantum vacuum, as expected. For r > 2 the line 
separating classical and quantum modes is approximately described by fc = *Jt. 



V. NON-LINEAR QUANTUM EVOLUTION AND SYMMETRY BREAKING 

To address the issue of symmetry breaking after inflation it is essential to incorporate the non-linear effects pro- 
portional to A. A full non-perturbative quantum treatment is beyond reach. However, we have seen in the previous 
section that the dynamics in the absence of non-linear terms gives rise to a fast growth of the amplitude of the 
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FIG. 2: The time for which a given mode k can be treated as classical |ffe(T c i)| = 1 is above the line in this figure. It is clear 
that long wavelength modes with < k ^ 1 become classical very early, at r c i — 2, while there remains, at any given time, a 
high energy spectrum of quantum modes, for k S> 1. 
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low- lying momenta, leading to wave- functions which are squeezed (quasi-classical in our language). We argue that 
even when the interaction is switched on the dynamics of these modes dominates the evolution of the system (at least 
during the first stages), and that this dynamics is described by classical field theory. The argument does not apply 
to higher momentum modes which sit largely in the quantum mechanical ground state. However, in quantum field 
theory, high momentum modes, although small, do not give negligible contribution to observables. Actually, naively 
their contribution is divergent. Nevertheless, we argue that the main contribution of the small quantum mechanical 
high-momentum modes sits in the renormalization of the constants to be used in the classical theory. 

It is possible to partially test this scheme in perturbation theory. Already at this stage the problem of infinities and 
renormalization arises In this section we will summarily analyse this issue, relegating the details of the calculations 
to the Appendix C. As we will see, for the program to be consistent one has to allow for a renormalization of the 
speed V of the inflaton at the bifurcation. 

In the standard setting, infinities in observables occur through the contribution of the infinite tower of momentum 
states. Introducing a cut-off in the problem makes the results finite, but cut-off dependent. It turns out, however, 
that in renormalizable theories, the only surviving effects of the cut-off at scales much smaller than itself are the 
modification of the constants of the theory. This allows the process of renormalization in which we recover uniqueness 
of the theory at the expense of taking this constants from experiment. We will now reexamine this problem for our 
time-de pen dent situation. Several research groups have investigated this problem in the past in different contexts, see 

Will. 

All the physical content of the theory is contained in the expectation values of products of the field operator at 
equal or different space-time points (we use the Heisenberg picture and expectation values should be understood as 
taken in the vacuum at r = tq): 

(y(n,xi) . . .y(r„,x„)) . (55) 

By differentiating with respect to r one can obtain expectation values of products of y and p. If we were to compute 
these quantities in the gaussian (non-interacting) theory, we would obtain, via Wick's theorem, a sum over all pairings 
of a product of factors associated to each pair 

G(°)(r,r',x-x')^( 2 /o(r,x) 2/0 (r',x'))= J ^e^'/.W/^r'), (56) 

where yo denotes the gaussian field for A = 0. The gaussian two-point function is the Fourier transform of £22 (fc, t, t') 
and is finite provided x/x' and/or r j^t' . 

The correlation functions at different times (i.e. Wightman functions) can be computed in perturbation theory 
by the method described in the Appendix C. Wightman functions are complex and unlike Feynman Green- functions 
(time-ordered products) depend on the order of the operators. According with our criterion for the gaussian case, 
we will consider Weyl-ordered (symmetrised) products to make the matching with the classical theory. If we now 
consider the symmetrised two-point function 

(y(x, r)y(x', r')) w = / e ik ( x " x > G(k, r, r 1 ) (57) 



(2tt) 3 

we can compute it to first order in A. The result is 

G(fc,r,r') = Re[/ fe (r)/ fe *(r')] + 2(N C + 2)A f ds A(s) Im [f k (r)f^s)} Re [/ fe (r')/ fc *(s)] (58) 



+ 2(N C + 2)X dsA(s)Im[f k (r')f^)]Mfk(r)m^] 







where iV c denotes the number of components of the Higgs field. The quantity A(s) gives the contribution of the 
tadpole sub-diagram, i.e. the two-point function at equal times and zero distance, and is given by 

A(t) = G(°>(t,t,0) = J ^l/ fe (r)| 2 = ^ / fp(k,r), (59) 

where the power spectrum is defined as P(k,r) = fc 3 |/fc(r)| 2 . This quantity is ultraviolet divergent. The structure of 
the divergence can be deduced by analysing the large k behaviour of the integrand. Using our previous expressions 
(with z — k 2 — t) and the asymptotic behaviour of Airy functions |3l| we get 



\Mr)\ 2 - ^ 



r / sin(2/cr) 



2k 2 V 2k 



T 



)+0{r 2 ) 



(60) 



P(k,r) = kV k (r)^ k 4 + Ul-^) + Or4). (61) 



2 4V 2kr 
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Thus A(t) has a time-independent quadratic divergence and a linear in time logarithmic divergence. 

Before explaining how can one deal with the divergence, we comment that Eqs. I|57I58|I coincide precisely with the 
calculation of the expectation values of the product of classical random field to the same order in perturbation theory. 
Divergences are hence present in both the quantum and the classical theory. Details of this calculation are also shown 
in Appendix C. 

We now address the problem of infinities that have occurred at this level. In the standard quantum theory the 
procedure is well-known. The calculation can be done using some regulator to cut off the contributions of high 
momenta, but this has to be accompanied by the addition of counter-terms in the interaction Hamiltonian. For 
the theory to be renormalizable these counter-terms should have the same expression as those appearing in the 
Hamiltonian (free or interacting) but with coefficients which are cut-off dependent and proportional to some power 
of A. This addition should get rid of infinities. Note that in our case a counter-term of the form 

_^+^ A ( 5l (A) + r5 2 (A))y 2 (x,r), (62) 

with <5i and 82 appropriately chosen cut-off dependent functions, is able to subtract the infinities encountered in A{t). 
Regularising the integrals by introducing a cut-off in momenta k < A we then get: 

A I < m (T,n)=A tes (T,A)-8i(A,( J ,)-T8 2 (A,(x) (63) 

To fix the arbitrariness introduced in the theory by the counter-term we must impose adequate renormalization 
conditions. As will be argued below, one convenient possibility is to choose the counter-term as 

8 1 (A,»)+t8 2 (A,v) = ^ J dkk 2 (\f k (T = fi 2 )\ 2 +(t-/! 2 )2F(t = m 2 )) . (64) 

where /1 denotes the characteristic mass scale of the problem which for a given time r is precisely y/r. We will call this 
renormalization prescription, the fixed-time subtraction scheme. Another possibility is a minimal subtraction scheme 
(not to be confused with the MS scheme of dimensional regularization) 

A rca (r, M) = A^(r, A) - JL (V + r log ±\ . (65) 

which differs from the previous one by finite terms of the form a + br. Actually, the renormalised quantity is obtained 
only after taking the limit A — * 00 in the subtracted quantity, but in practice taking A sufficiently large is a good 
approximation. 

The fact that the structure of the counter-terms (or of the divergence) has the same form as the terms already 
present in the Hamiltonian, shows that our calculation is consistent at least to this order. The infinities are re-absorbed 
in the values of the constants of the theory. Here, in addition to the ordinary time-independent subtraction we have 
a counter-term linear in r, which can be interpreted as a renormalization of x — Xc, i-e. of the inflaton velocity V. A 
different choice of scheme is compensated by a finite renormalization of the parameters of our model. 

Now we look back at the problem of approximating the result by a classical random field. Since the regularised 
result to this order is the same (for symmetric expectation values) a similar subtraction procedure is necessary. There 
is certainly no problem to do so in perturbation theory. However, in practice what we want to do is to be able to match 
the renormalised quantum result by modifying the initial spectrum of the classical field to be used as starting point 
for the classical evolution. Notice that when r = Tj = /1 2 the value of A IQn obtained with the fixed-time subtraction 
scheme, 11640 . is exactly reproduced by truncating the initial spectrum at /1 = yfri. This is a very natural choice from 
the point of view of the classical approximation. As can be seen from Figs. ^ 121 f°r large enough r the separation 
between quantum and classical modes sits indeed at k ~ \fr. At a given time modes with momenta below yfr have 
been amplified while those above \fr remain in the vacuum. The amplification proceeds until some time r s b, when 
(0 2 (r s b)) gets close to the vacuum expectation value v 2 and the field starts oscillating around the true vacuum. The 
dynamics of symmetry breaking is hence expected to be governed by the low momentum modes with k 2 < r s b whose 
evolution can be described in the classical approximation (as we will see below for a large range of parameters, r s b 
varies only within the values r s b = 5 ± 2). The classical theory can then be seen, in a way analogous to what happens 
at high temperature 35, 36], as an effective theory where momenta above £;* = ^/r s b have been integrated out. As 
far as modes above fc* are not highly populated by rescattering and back-reaction this effective theory is expected to 
be valid and can be studied within the classical approximation. 

In summary, our proposal is to fix our classical field by matching its correlation functions with the renormalised 
perturbative expression at a time Tj = /1 2 such that a sufficiently large number of momentum modes have become 
classical but well before non linearities have set in. The initial spectrum of the classical field will be cut-off at 
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FIG. 3: The power spectrum of the Higgs quantum fluctuations, P(k,r)/k 2 = k\fk{r)\ 2 , at different times in the evolution. 
The dotted vertical lines indicate the value of the cut-off, at k — \pr, where the classical spectrum is truncated. Also shown is 
the excellent approximation 1691 in the region of long wave modes. 



k = k„ = ^fri. This eliminates the UV infinities of the classical theory. If we compare now with the calculation at one 
loop, we realise that the parameters entering the classical theory are the renormalised parameters in the fixed-time 
subtraction scheme <|64|) . As we will see in what follows and in Section lYl Bl our results are fairly insensitive to the 
specific choice of Tj within a scaling window below r s b- 

The validity of this approximation can be partially tested in perturbation theory. A first check is the form of 
the power spectrum (|59|) . We have plotted in Fig. [21 the power spectra P(k,r) divided by k 2 , for four values of the 
normalised time r — 2, 3, 4, 5. Clearly, the power spectra grow in time faster than an exponential, at a very large rate 
in fact. We take as initial spectrum of the classical field at a given initial time Tj, the exact power spectrum cutoff at 
k = \fri- As seen in the figure this encompasses almost all the physically relevant low momentum modes for T{ ^ 2. 

In Fig. 21 we compare, for several values of Tj, A IBa (r, \i = y/fi) in the fixed-time scheme with A claB (r, r^), obtained 
from cutting off the power spectrum, P(fc,Tj) in Fig. |3 at k = <Jri. |40) . For Tj = 2 the maximal difference between 
them amounts to 2%, rapidly decreasing as we increase Ti. A direct comparison between the values of A r<m and A claB 
for Ti — 2 is also shown. The goodness of the approximation performed by truncating the spectrum is clearly evident. 
We also study the dependence of A rea on the value of Tj used for the fixed-time renormalization scheme. We plot the 
difference between A IO „ defined at = 2 and Ti = 3. As it should, it is of the form a + br and it remains very small 
in all the range of times we are interested in. 

It is easy to estimate the time r s b when symmetry breaking is expected and the amplification of modes ceases to 
take place. We can estimate the time of symmetry breaking r s b by equating 

(|0| 2 (r sb )) = M 2 N c A tm (T sh , i± = y/ri) = v 2 = N c v 2 . (66) 

We have just described how for Tj 2 a very good approximation for A ron , in the fixed-time subtraction scheme, is 
obtained by just truncating the power spectrum at /x 2 = t^. We can thus approximate the above expression for the 
vev by 

M 2 N 

<H 2 )-^/ o dke\f k \ 2 = v 2 (67) 
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and fSty . 



which can be rewritten as 

rs/n (ifo 2tt 2 2tt 2 

P{Tsh) "Jo T P{k ' Tsh) = XNjwW S A WW 3 (68) 

This can be computed exactly using (|43|) but to give an analytic estimate of its dependence on the parameters, we 
will use an approximation to the classical power spectrum (|59|l . In the region of quasi-classical modes it is very well 
described by, 

P apv {k 1 r)=A{r)k 2 e- B ^ k \ (69) 
A{t) = A Bi 2 (r), B(T)=2Vr-2, (70) 

which can be obtained from \fk\ 2 — |Ci(fc)| 2 |Bi| 2 (z), where 

fc|Cl(fc)|2 "Tr^w( 1 + 2fc2 + 0(fc4) ) " Aoe2k2 > (71) 

exp (Jz 3 / 2 ) = exp (Jr 3 / 2 - 2^/7 fc 2 + C(fc 4 )) . (72) 

We have plotted P app (fc, t) together with the exact spectrum in Fig. 01 We can see that it is an excellent approximation 
to the classical power spectrum, in the region of interest. Using P app (fc,-r) to estimate (16811 gives the condition 

( \ 2n2 f dk p ru \ ^(^b) 

p(Tsb) = A (2V)V3 J T Fapp(fc ' Tsb) = 25M ' (73) 
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TABLE I: The time scales of symmetry breaking and the onset of the non-linear stage for different model parameters. The 
coupling depends on the number of components of the Higgs field as Aq = N C X. 



V 


Ao 




mi n i 


Tab 


mt sb 


0.003 


0.11 


2.78 


15.3 


4.6 


25.3 


0.003 


0.01 


3.76 


20.7 


5.2 


28.5 


0.0003 


0.001 


4.01 


47.5 


6.0 


71.3 


0.00002 


0.0001 


4.82 


141. 


6.8 


200. 



We have evaluated this function p(r s b) numerically and found an excellent fit to it, in the range r > 1, as, for N c = 4, 

lnp(r s b) = -3.5+[8 + r s 3 b 23 ] 1 /2, (74 ) 

which gives directly the time of symmetry breaking in units of m , 

mtsb = (2^)- 1/3 

We can use this compact expression to estimate the time of symmetry breaking for any coupling A and any inflaton 
velocity V at the bifurcation. For example, for Ao = 0.11 and V = 0.003, we find r s b = 4.6 and mi s b = 25.3, which 
agrees very well with numerical (lattice) simulations performed for those values of the parameters. 

Note that, as mentioned before, the dependence of r sb with the parameters A = \q/N c and V is very mild. In the 
whole range of parameters, Ao (2V) 2 / 3 € [10 -8 , 1], the normalised time of symmetry breaking only varies within the 
range r s b = 5 ± 2. Some particular examples can be found in Table [I] 

VI. NON-LINEAR EVOLUTION OF THE CLASSICAL SYSTEM 

In the previous sections we have argued that the first stages of the quantum evolution of the system (when the 
non-linear self-coupling of the Higgs is negligible) drive the system into a state with highly populated low momentum 
modes. The evolution of this state can be accounted for by the evolution of a classical (approximately gaussian) 
random field. This justifies the main assumption of this and remaining sections, namely that the subsequent non- 
linear dynamics of the system is determined by the classical evolution of this field. This evolution is deterministic and 
the random character appears in the initial values of the field at time r = Tj. These initial conditions are determined 
by the exact gaussian quantum evolution of the system studied in the previous section. Thus, the initial Higgs field 
is chosen gaussian, an approximation which can be tested by probing the sensitivity of our results to the value of Tj. 
As we will see this works very well within the appropriate range of initial times. Some statistical properties of this 
initial gaussian random field can be studied analytically. This is done in section 6.1. These properties extend to times 
during which the evolution is essentially linear and the field remains approximately gaussian. A full non-perturbative 
treatment of the dynamics can only be done by numerical methods. We have actually carried this step by lattice 
simulations. This is described in section 6.2 where a full account of the methodology and the checks performed to 
show cut-off independence is described. Results will be presented in the next section 

A. Peaks of the Higgs spatial distribution 

The statistics of the Higgs spatial distribution can be determined from the gaussian fluctuations that are used to 
build it up. A detailed description can be found in Ref. [24| for the case of the gaussian density field responsible for 
galaxy formation. In fact, the spatial distribution and subsequent dynamical behaviour of the Higgs field at the initial 
stages of symmetry breaking turns out to be not that different from that of both the linear and non-linear growth of 
the cosmological density field (also built up from the gaussian random fields of cosmological perturbations), except in 
the dynamics of gravitational collapse of the latter. 

The fact that the quantum fluctuations of the Higgs give rise to a classical gaussian random field, allows us to study 
the statistical properties of this field in terms of a single function, the two-point correlation function in Fourier space 
(i.e. the power spectrum), whose approximate expression can be found in (|69|) . This quantity allows the computation 
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of several related quantities that characterise a gaussian random field, e.g. the spatial correlation function, the density 
of peaks above a certain threshold, the shape of the highest peaks, etc. 

The first quantity that we can compute is the spatial correlation function, defined as the two-point correlation 
function between two points separated by a distance r, 

„ s , ,/ s ^ M 2 N C f°° dk „ . sin kr 



2tt 2 L k " rov '"' fer 



M 2 N c A(t) 
2tt 2 ~ 
Af 2 7V c ^(t)^ 



dke- B{r)k2 sin fc 



r 



where erfi(a;) is the imaginary error function j31j . This correlation function determines the average size of the lumps, 

6>, 



Zo~2B 1 I 2 {t)~2^2^-2, (77) 

as a function of normalised time r. Note that the time dependence of the correlation length is different than for a 
quench symmetry breaking. While in the latter case, the correlation length grows like £o ~ 2-y/r, in our case, it grows 
like £o ~ 2V2T 1 / 4 for "large" r (still in the linear regime). This introduces some slight differences in the behaviour 
of the field at symmetry breaking. 

We can then compute from (|76l) the time-dependent dispersion 

a(r) e, e /2 (0,r) = rms (r) = ^^ p V\ T ) , (78) 

V I 7T 

which is nothing but the root mean square value of the Higgs field. 

Another quantity which is very useful to characterise the field distribution is the number density of peaks of the 
field above a certain threshold <j) c , see Ref. pil W\ 



WW = — 2 ) {V - 1} 6XP( ~" /2) ' (79) 

-£"(0,r) _ (k 2 ) J(dk/k)P(k,r)k 2 



npeak(r) = (y 2 - 1) exp(-^ 2 /2) , (81) 



£(0,t) 3 3j(dk/k)P(k,r) ' 

where v = > 1. In our case, the number density of high peaks is given by 

1 

3^3 tt 2 rf(r) 

( \ ^ c ^ 2_7r -1/2/ \ /o \ 

v[t) = ■ — 1 (r . (82) 

V ' v (2F) 1 /3yA K ' ^ ' 

We have evaluated this function for the parameters A = 0.11/4 and V = 0.003, at various times and compared with our 
lattice simulations (for different volumes V — (27r/p m i n ) 3 and lattice spacings a). The results are very encouraging. 
If we multiply this density of peaks by the actual volume of the simulations, we find indeed just a few peaks above 
e.g. 4> c — 0.02w, at the time of symmetry breaking. 

In fact, we can compute not only the probability per unit volume to find a peak in the distribution of the Higgs 
field, but also their radial profile |24j . 

, , MJW C f°° dk 1/9/ , . sin kr 
p(r,r) = / —pV 2 (k, T ) 



V2ir Jo k rcg V ' ' kr 
M^W c A^ 2 (t) f°° dk 



V2nr Jq k 



e 2 sin kr 



2V2r \^2B(t)) 
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FIG. 5: The radial profile of the Higgs peak for A — 0.11/4 and V = 0.003, at time r = 2.54, corresponding to rat — 14, obtained 
with our lattice simulation (with error bars, from averaging over several realizations), and compared with the analytical result 
1881 . We have also included the rms Higgs value j^Bj a t that time. Note that we are still in the linear regime, where l|83|l gives 
a very good approximation. The higher tail corresponds to an averaging out of several lower peaks. 



where erf (2) is the error function [31]]. We have plotted this profile function in terms of the radial coordinate, together 
with the lattice results in Fig. for A = 0.11/4 and V = 0.003, at time r = 2.54, corresponding to mt — 14, well 
before symmetry breaking, which occurs at mi s b — 26. 



B. Lattice simulations 



The previous analysis falls short of addressing the most important aspects of symmetry breaking after hybrid 
inflation since the main effect is non-perturbative, see Ref. |lj|- As discussed we will incorporate these effects by 
performing classical real-time numerical simulations in the lattice. Generically this classical approximation would 
fail to reproduce the relevant physics but we have just argued that this is indeed the correct approximation for the 
infrared modes of the Higgs at the time of symmetry breaking. 

The usual procedure |38| |41| is to take as initial conditions for the lattice simulations gaussian random fields given 
by the distribution l|34|) with vacuum initial amplitudes corresponding to Eqs. I|17|) . (|18fl . We would like to stress 
here that the correct description of the quantum linear system in terms of a gaussian random field requires the use 
of two independent gaussian variables, as indicated in l|34|l . One of them, y in 1341) . describes field fluctuations with 
dispersion \fk\ 2 and a random phase. The other, p, with dispersion (4 |/fe(r)| 2 ) -1 and a random phase, allows to 
define the conjugate momentum through, see Appendix B, 

l/feO)r 

Notice that this prescription is valid in order to give initial conditions at any time during the evolution before 
non-linearities set in. In particular, as described in section we propose to take as starting point for the lattice 
simulations the above gaussian ensemble at a fixed time sufficiently advanced to guarantee that a large fraction of 
modes have become classical, but well before the time when non linearities become relevant, in a similar spirit as that 
in Ref. Q. This has the advantage of allowing a clear separation between infrared (classical) modes which evolve 
classically and ultraviolet (quantum ) modes that will be absorbed in the renormalization of the constants of the 
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theory. From the previous analysis, see Fig. [2 and the discussion after Eq. Q75[l. a good choice for the matching time 
in a wide range of model parameters seems to be r, ~ 2. See the discussion in the next section about the onset of the 
non-linear regime. 

Therefore, we propose the following as initial conditions in our lattice simulations. At a fixed time Ti previous to 
symmetry breaking: 

• Put to zero all the modes that have not become classical at r^. This includes all the modes of the inflaton but 
the homogeneous zero mode, and all large momentum modes of the Higgs with k > ^fri (replacing the hard 
cutoff at yfri by the approximate power spectrum in Hfi9[). which strongly dumps ultraviolet modes, does not 
significantly change the results even at a quantitative level). 

• Set the homogeneous zero mode of the inflaton to x/Xc = 1 — Vmti with conjugate momentum x/Xc = —Vm. 

• For the Higgs fluctuations, each Fourier component, with momentum \k\ < y/ri, has an amplitude \<pk\ randomly 
generated according to the Raylcigh distribution: 



with dispersion given by a\ = \ fk\ 2 = k P(k,Ti), and a uniform random phase 9k € [0, 2w]. The conjugate 
momentum gk = i(j)' k is uniquely determined once <f>k is known, through the relation 



with fk(n) and Fk(Ti) given by Eqs. (|4*5|) and Q26|). respectively, at r = 7j. This corresponds to the classical 
limit of (|84|) . an approximation that is well justified for ^ 2, see Appendix B. 

• Take the masses and couplings used in the simulation as the physical renormalised ones in the fixed-time 
subtraction prescription. 

As long as the time chosen for initialisation is sufficiently advanced that a large fraction of modes have become 
classical, we hope that most of the physics responsible for symmetry breaking will be included in the simulations. 
How advanced it has to be in a concrete realization can be tested by studying in which range the time evolution is 
insensitive to the choice of r^. This provides also a check of the validity of our approach. The result of such a test 
is presented in Fig. We compare the time evolution of (\4>\}, (x) & n d the average kinetic, gradient and potential 
energies obtained from setting the initial conditions at mti = 10,12,14 (r^ = 1.81,2.18,2.54) for A = 0.11/4 and 
V = 0.003. The agreement is excellent, corroborating our estimate that for Ti ~ 2 all the basic relevant modes driving 
symmetry breaking have already become classical and thereafter the evolution is well described by our lattice classical 
simulations. 

All the lattice results presented in this paper have been obtained for a SU(2) Higgs doublet coupled to the inflaton 
with coupling g 2 = 2A = 0.22/4, and inflaton velocity V = 0.003. Due to the finite volume, the momentum k 
is discretised in units of a minimal momentum given by p m i n = 2~k/L, with L = Na, where N is the number of 
lattice points, and a the lattice spacing. Our simulations have been performed in lattices of sizes 32 3 , 48 3 and 64 3 
with physical volumes determined by p m i n — 0.1m, 0.075m and 0.05m and lattice spacings varying from ma ~ 1 to 
ma ~ 2. The choice of lattice volumes and lattice spacings has been performed such as to avoid lattice spacing and 
finite volume dependence of the lattice results. Notice that the minimal momentum has to be small enough that a 
sufficiently large number of classical momenta with k 1 is taken into account. We have found that for p m - la < 0.1m 
this is indeed the case and no significant volume dependence is observed. 

A further essential test of our approach is that it succeeds in taming ultraviolet divergences. On the lattice there 
is a maximal momentum determined by the lattice cutoff through p max = 2ix j a. Naturally, re-scattering and back- 
reaction will populate the high momentum modes at and after symmetry breaking This is certainly a physical effect 
but if the lattice cutoff is not chosen large enough population of the high momentum modes is artificially induced 
by cutoff effects. The lattice cutoff should then be chosen such as to avoid that this takes place before the relevant 
dynamics of symmetry breaking. A reasonable value for our choice of parameters is ma < 2, as can be seen from 
Fig. El where we compare the time evolution of (\(j>\), (x) an d average energies for several values of the lattice cutoff 
ma = 0.98, 1.31, 1.96. No significant lattice spacing dependence is observed here, while it becomes clearly appreciable 
for ma > 3. Details of the simulations and further results will be presented elsewhere. These lattice simulations will 
allow us to test the next stage, the non-linear approach to symmetry breaking. 
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FIG. 6: The time evolution of (\<j)\), (x) an< l energies (normalised to the initial one) for A = 0.11/4 and V = 0.003, obtained with 
our lattice simulation. Left: for different choices of mti, the time for matching the quantum evolution to the classical lattice 
simulations. Right: for different values of the lattice spacing ma — 0.98, 1.31, 1.96. 



VII. RESULTS OF THE NON-LINEAR ANALYSIS AND "BUBBLE" FORMATION 



In the previous sections we have developed a formalism to describe the linear growth of the Higgs quantum fluc- 
tuations and their conversion into a classical gaussian random field. As we have argued, in the linear regime there is 
only the homogeneous mode of the inflaton, (x/Xc) — 1 — Vm(t — t c ), which induces a negative mass squared |J5J for 
the Higgs, and thus its spinodal instability towards the true vacuum. The quick growth of the quantum fluctuations 
generates a gaussian random field with correlation function (|7fj|) . and a rms field value, see (J7SJ, 



rms (2I/) 1 / 3 VAlVe 1/2 , , 

(v = — = /= P ' 0) • 

V V2 7T 



(87) 



where p{r) is given by (|68|l . Eventually the mean field <j> will become large and will approach the vev of its potential, 
thus breaking the symmetry. Before that happens, its coupling to the inflaton will induce a back-reaction on the 
homogeneous mode of the inflaton, \, which will start to deviate from the linear regime described above. At this 
stage the non-perturbative evolution can be studied by numerically solving the coupled classical equations of motion 
for the inflaton and Higgs: 

a (x, t) - V V(x, t) + (|0| 2 (x, t) + X 2 (x, t) - 1) a (x, t) = , 

2 

X (x, t) - V 2 X (x, t) + ^- H 2 (x, t) X (x, t) = , 



(89) 



with gaussian initial conditions as described in the previous sections. In this section we will present the results of our 
numerical simulations and give an approximate analytic understanding of how symmetry breaking takes place. 
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mt = 23 mt = 24 




FIG. 7: Snapshots of the growth of the Higgs peak in a full non-linear lattice simulation, for A = 0.11/4 and V = 0.003. Plotted 
is the value of the Higgs amplitude (j> in the plane (x, y), where the ^-coordinate is that of the highest peak. Note that several 
peaks appear in the simulation. Here we show the first stages of the evolution, where the highest peak invaginates and forms 
what we call the "bubble". 

Although the initial conditions are random, as a result of the non-linear dynamics many of the qualitative features 
of the evolution are fairly universal, although quantitatively different configurations differ by small shifts in the origin 
of times as well as spatially random positions for the center of the peaks. Therefore, we prefer to illustrate our analytic 
formulas by comparing with the results of a typical lattice configuration, e.g. the one displayed in Figs. 16181 

Symmetry breaking in our model is not at all a homogeneous process. Already in the linear regime, the Higgs 
field evolves by developing lumps in space that grow with time, see (|83f) and Fig. [5] The classical evolution of the 
Higgs' lumps, once non-linearities become relevant, can be followed in Figures {7\ and [S] where we show some snapshots 
of the growth of the Higgs' peaks from the first stages of the evolution, rat = 23, till mt = 40 above which full 
symmetry breaking takes place and the mean Higgs field approaches the vev. As can be seen from the figures, the 
peak of the largest Higgs' lump is the first to break the symmetry, i.e. to reach \(f>\ — v, and soon after the center 
of the lump invaginates, creating an approximately spherically symmetric "bubble" , with "ridges" that remain above 
1 4>\ = v. Finally, neighbouring bubbles collide and the symmetry gets fully broken through the generation of higher 
momentum modes. In Fig. [!|]we show the behaviour of |^>(x, t)\ at the center of the highest Higgs lump. It oscillates 
around \4>\ = v with an amplitude that is dumped in time. Oscillations remain coherent giving rise to concentric 
bubbles, until the time when bubble collisions break the symmetry. 

It is possible to get an analytic understanding on how this non-linear process takes place before bubbles start to 




FIG. 8: Same as in Fig.Q Here we show the late stages, in which gradients arise from collisions of bubbles and the symmetry 
is broken, i.e. <j> — 1- 



collide. For the problem we are considering, we can rewrite the components of the Higgs field as (jf = <fi h a (we will 
use from now on the symbol <fi to denote the modulus of the Higgs) while ft = h - <r £ SU(2) is an element of the gauge 
group, with a = (11, it) with r a the Pauli matrices. With this the equations of motion for the coupled inflaton Higgs 
field can be rewritten as 

0(x, t) - V 2 0(x, t) + {4> 2 + x 2 - 1) 4> - l^n| 2 - , (90) 
d' i {(j> 2 d ll n)=-n(j> 2 \d lx n\ 2 , (91) 

X(x,i)-V 2 x(x,t) + ^ 2 X = 0, (92) 

where dots and V denote derivatives with respect to mt and mx respectively, and the homogeneous modes have been 
normalised to their vevs, <fi/v — » <j> and x/Xc — * X- 

We can take advantage of the fact that, for g 2 = 2A, a solution to the set of coupled equations of motion is given 
by <f>(x.,t) = 1 — x(x, t) \l4, |26I| and d^{4> 2 d^n) = 0. Numerical results corroborate that this is very approximately 
the solution soon after non-linearities set in. In Figure we show, for our model with parameters V — 0.003 and 
A = 0.11/4, <fi versus 1 — x at the location of the highest Higgs lump. Comparing with Figure |H1 we can follow how 
the Higgs and inflaton evolve colinearly during all the time of coherent oscillations of the peak. In Fig. 1101 we show 
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FIG. 9: Left: The time evolution of the modulus of the Higgs, at the location of the highest Higgs peak. Plotted is the Higgs 
modulus 4>/v as a function of time. Note the effect of bubble collisions on the Higgs oscillations after mt = 40. Right: Collinear 
evolution of the inflaton and the Higgs at the location of the highest Higgs peak. Note that the inflaton and Higgs satisfy 
<j> = 1 — x to very good accuracy, until rather late, when "bubbles" start to collide. 
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FIG. 10: The histogram of field values for the Higgs modulus <f> and the co-inflaton 1 — \. Note that soon after the initial 
condition and the subsequent non-linear evolution, the two coincide very precisely. Note that a few oscillations can be observed 
during symmetry breaking, i.e. during times mt — 22 — 30. 



the distribution of field values as a function of time, in the time interval between mt = 12 and mt = 60, where most 
of the action takes place. During most of the non-linear initial stage, through symmetry breaking and until bubbles 
collide we have: 0(x, t) = 1 — x(x, t), Vx . 

During the time that inflaton and Higgs evolve colinearly, the system can be seen as that of a single field with a 
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FIG. 11: The time evolution of the Higgs peak (r = 0) and the Higgs rms value, obtained with our lattice simulations, as 
compared with the analytical result (19711 . and the numerical solution of equation 11001 . which includes the gradient terms. Also 
shown is the comparison between the Higgs and the inflaton evolution, i.e. <f>(t) and 1 — x(t)- 



modified potential V (</)), with the minimum at cf> = 1. 



1 



-(d^) 2 -v{<t>) 



V{4>) = ^{1-^ + 3^ 

The equation of motion of the scalar field <f> becomes a non-linear partial differential equation 

4>(x, t) - V 2 c/)(x, t) - 2c/) 2 (x, t) + 2c/> 3 (x, t)=0. 



(93) 



(94) 



If the gradient terms are much smaller than the non-linear ones, we can as a first approximation neglect them leading 
to 



0(f) - 2c/> 2 (l - (f>) =0, 
which leads to a conserved energy E = Eq + | with 



(95) 



(96) 



A solution with Bo = 0, a very good approximation taking into account that initially both the field and its derivative 
are very small, is given by 



<t>(t) = 1 - x (t) 



12 



9 + 4[mi-mi max ] 2 



(97) 



with i max the time at which the field reaches its maximum value: c/)(i max ) = |. This time can be rewritten in terms 
of the value of the field at any other time t n \ through 



3 9 

h ~ 4 
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FIG. 12: The time evolution ol the Higgs radial profile around the highest peak, obtained with our lattice simulation (points), 
as compared with the numerical solution of the partial differential equation (lines). It is surprising how well the formation 
of the bubble is reproduced with the simple assumption of homogeneity around the peak of the bubble. Of course, the peak 
solution does not take into account the presence of secondary bubbles, that appear in the lattice simulation at r ~ 30. 



In particular, at every point x, we can take <po as the "initial" value of the Higgs field. This is given by the profile 
of the lump in the linear approximation, Eq. (|83l) . at a time t n \ at which the evolution becomes non- linear and we 
can no longer ignore its higher order interactions. In Fig. 1111 we show again the non- linear growth of the Higgs field 
at the top of the largest peak in the simulation, and compare it with the analytical solution (|97|) . The agreement 
is very good during the first oscillation although (|97|) cannot reproduce the subsequent ones. At these stage we can 
already understand how the spherical bubbles arise. Take the spherically symmetric peak profile (|83() at the non-linear 
time t n i and let each point x evolve like l|97[l . Points with higher value of 4>o(t n i) will reach first the maximum value 
(0max = § ) and then decrease. This generates a spherical wave that propagates from the center of the lump to infinity. 
The production of bubbles associated with symmetry breaking were first described in Ref. [l4T | for the model A</> 3 , 
which is analogous to our reduced model. The subsequent evolution is of course different, due to the presence of the 
infiaton field. 

We can evaluate the non-linear time, t n i, by equating 



X (r a {) = l-\{2V) 2 l\ nX 



l - Mtox) 



Using H87f) and (|74ll . we can find the nonlinear time t„i as the solution of the transcendental equation 

, 0.31 



r 



3. 5 + 2 



(99) 



For the values of parameters chosen, A = 0.11/4 and V = 0.003, we find mt n \ = 15.3 and 0o = 0.1. That is, soon 
after the Higgs field becomes non-linear, it ceases to grow exponentially like Ij87(l . and starts to grow like (|97|l . which 
has a peak at mt s ~ 23 J$ mtsb ~ 26, see Fig. 1111 This corresponds to a time slightly earlier that the time of 
symmetry breaking. This is of course natural since, as we have described, the Higgs field has an inhomogeneous 
spatial distribution. The mean field (coarse-grained over a horizon-sized volume) is much lower than a typical peak 
of the field. The top of the peak follows very approximately the homogeneous equation (|95|l , with solution (|97|) . High 
peaks will reach the symmetry breaking vev much earlier than the mean field, and will oscillate around the vev with 
a much larger amplitude that the average (coarse-grained) field. 
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Obviously, the phenomenological damping of oscillations that we have described has to arise from the gradient terms 
which we have neglected. Hence, we will improve our approximation by keeping these terms, but assuming spherical 
symmetry (<jf>(x, t) — > (j)(r, t)) around the center of the lump (r = 0). Our lattice data support the approximate validity 
of this assumption. This will allow us to track the time-evolution of the lump profile as it develops into bubbles. The 
two-dimensional partial differential equation for 0(r, t) becomes: 

6(r, t) - 4>"{r, t) - -6'(r, t) - 20 2 (r, t) + 20 3 (r, t)=0. (100) 
r 

We have solved this equation numerically. The initial condition was fixed at a time mt n \ when the profile matches 
expression |J53). In order to compare with the non-linear lattice simulations, we added by hand a tail a long distances 
, to match the lattice initial conditions, see Fig. [3] To fix a unique solution, one has also to fix <j)(r, t n \)- Choosing 
this derivative equal to zero we have obtained the shape 4*(r,t) for all times in the region of interest. In Fig. 1121 we 
present the result of the comparison of these results with those obtained from the full 4D lattice real-time equations 
of motion. The general shape is quite properly reproduced. Furthermore, the oscillations of the peak height are also 
recovered, see Fig. \H\ 

As a last remark, note that the bubbles that appear here are not vacuum bubbles like those produced in a first order 
phase transition, since the interior of them is not in the true vacuum. Furthermore, we note also that the ridges of the 
bubbles are moving very fast and presumably subsequent collisions between bubbles formed at different space-time 
points arc highly relativistic, and may be responsible for a large density of gravitational waves, which could be seen 
in LISA. 

For a typical lattice configuration one can follow the evolution of the Higgs from the formation of the 
first bubbles to the breaking of the symmetry with the .gif file that can be found in the web page: 
http://lattice.ft . uam . es/SymBrk/2dHiggs . gif 

VIII. CONCLUSIONS 

In this paper we have studied the evolution of a hybrid inflation model from the quantum false vacuum state at 
the end of inflation to the broken symmetry true vacuum state. A full description of this dynamics amounts to a 
non-perturbative, non-linear, real-time evolution of the quantum system, which looks a priori like a formidable task. 
The size of non- linear effects is given by A0 2 , where A is the coupling constant and <f> 2 the square of the typical value 
of the Higgs field. Since initially <f> 2 and A are small, it is reasonable to assume that perturbation theory is a good 
approximation and the dynamics is well approximated by the gaussian Hamiltonian. However, the quantum evolution 
of this gaussian system, which can be treated exactly, is far from trivial. The complexity results from the negative 
time-dependent mass-square of low-momentum Higgs modes induced by the coupling to the inflaton. This tachyonic 
dynamics generates a faster than exponential tachyonic growth of low- lying momentum modes of the Higgs, giving 
rise to regions where A</> 2 is non-negligible and where non-linearities set in. In this paper we have shown that the 
dynamics of the tachyonic modes is well-described by that of a classical gaussian random field, a result that holds even 
after including perturbative corrections in the coupling, which are still accessible to exact computation. At this stage 
important considerations set in through the appearance of ultraviolet divergences. High-momentum modes cannot be 
neglected but their effect can be absorbed in the value of the couplings of the theory. Here, in addition to the usual 
standard time-independent renormalization, a renormalization of the initial velocity of the inflaton field is required 
to get rid of the time-dependent infinities generated at first order in the coupling A. 

The previous analysis justifies the next stage of the study carried out in this paper, namely the classical non-linear 
evolution of the resulting classical field. This problem can be addressed numerically by formulating the problem on a 
spatial lattice and evolving the system according to the classical real-time evolution equations. The initial conditions 
on the classical field are determined by the previously computed (non-self-interacting) quantum Higgs evolution. Our 
results are independent of all cut-offs introduced by this numerical procedure: the initial time of the simulation, the 
lattice-spacing and the finite lattice volume. This, of course, provided they are taken in the appropriate ranges. 

The resulting non-linear evolution which drives the system towards symmetry breaking is fairly non-trivial. The 
inhomogcncous Higgs field distribution has lumps in space whose height grows with time during the approximately 
linear evolution phase. This growth continues, although at a slower pace, when the non-linear terms become relevant. 
The behaviour changes again as the highest lumps reach the magnitude of the Higgs vacuum expectation value. 
Then the lumps evolve into approximately spherically symmetric bubbles which expand at a very high speed. It is 
important not to confuse these bubbles with those appearing in a first order phase transition which separate two 
different phases. Our bubbles are rather like spherical shock waves as those appearing in Ref. [1J| . These stages of the 
non-linear evolution can be qualitatively and quantitatively understood analytically. The last phase of evolution arises 
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as neighbouring bubbles collide and generate higher momentum modes. This phase is harder to tackle analytically 
but its early stages, at least, seems relatively safe for our lattice numerical procedure. 

In the early stages of evolution our results resemble those obtained for a one component Higgs model in Ref. (lfij . 
as expected from the decoupling of the different components of the Higgs field in the linear regime. At later times, 
however, the comparison is difficult due to the different nature of the defects in both theories. 

The authors are presently studying how the previously described processes might be influenced by the coupling to 
gauge fields, and its application to the study of physical phenomena such as baryogenesis. We anticipate that there is 
no essential obstruction for incorporating gauge fields, although the formalism complicates considerably. Furthermore, 
the numerical evolution including gauge fields does not change substantially the gross features of the picture described 
here. All this will be the subject of a future publication. 
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APPENDIX A: THE FORMALISM OF SQUEEZED STATES 



In this appendix we will summarise the concept of squeezed states so often used in quantum optics, and recently 
applied to the study of quantum fluctuations from inflation 0, 0, [2(j . 

The canonical harmonic oscillator system is described by two complex functions (/&, g k ), plus a Wronskian 
constraint l|19l) , and thus we can describe the system in terms of three real functions in the standard parametrisation 
for squeezed states, 

u k (r) = ^=[fc/ fc (r)+ 5fc (r)] -e-' i9fe(T) cosh r fe (r), (Al) 
V2fc 

v k (r) = -L [k ft (r) - g* k (r)] = e 4 e « ^ +2t « sinh r k (r) , ( A2) 

V Ik 

where r k is the squeezing parameter, (j) k the squeezing angle, and k the phase. 
We can also write its relation to the usual Bogoliubov coefficients, {a k , f3 k } 

u k = a k e~ lfcT , v* k = fa e lfcT , (A3) 

which is useful for the adiabatic expansion, and allows one to write the average number of particles and other quantities, 

n k = |/? fc | 2 = M 2 = 1 | 5fe | 2 + k - |A| 2 - 1 = sinh 2 r k , (A4) 

a k = 2Re (a* k fa e 2lfer ) = 2Re (u* k vfj = cos20 fe sinh2r fc , (A5) 

r k = 21m (a* k f3 k e 2ifcT ) = 21m (u% v£\ = - sin 2cj) k sinh 2r k . (A6) 
We can invert these expressions to give (r k , 9 k , cf> k ) as a function of u k and v k , 



sinh r k = y Re v k + Im v k , cosh r k — y Re u k + Im u k , (A7) 

tan 9 k = - ^ , tan(0 fe + 2<f> k ) = ^ , (A8) 

Re u k Re v k 

Im v k Re u k + Im u k Re v k 

tan20 fe = . A9 

Re v k Re u k — Im u k Im v k 

Let us now use the squeezing formalism to describe the evolution of the wave function. The equations of motion 
for the squeezing parameters follow from those of the field and momentum modes, 

w' 

r' k = — cos20 fc , (A10) 
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w' 



k coth2r^, sin20j. , (All) 



w 
w' 

6' k = k H tanh 2r fe sin 20 fc , (A12) 

w 

where we have replaced the time-dependent mass © with the function w, with 

w" =tw, — ► w{t) = Bi(r) + V3 Ai(r) , (A13) 

with Ai and Bi the two independent Airy functions [3l|], satisfying w'(0) = 0. 

As we can see in Fig. 1131 the evolution is driven towards large ^> 1. Thus, in that limit, 

la w ' sin 2 ^fc n 

(Vk + (Pk) = r-rr, > 0, 

w sinn Irk 

and therefore Or + 4>k — * const. We can always choose this constant to be zero, so that the real and imaginary 
components of the field and momentum modes become 

fki = -j= e rk cos0 fc , f k2 = -j= e~ rk sin0 fe , (A14) 



9ki = y ^ e rfc cos0 fc , 3fc 2 = y ^ e rfe sin0 fc . (A15) 

It is clear that, in the limit of large squeezing (r^ — > oo), the field mode fk becomes purely real, while the momentum 
mode gk becomes pure imaginary. This means that the field and momentum operators (|2U[1 become, in that limit, 



y(Kr) - V2kMr)y(Kr ) 1 ^ ^ f^-y^r). (A16) 

p(k,r) -» V2fc 5 fc2(r)y(k,r ) J /*i( r ) 
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As a consequence of this squeezing, information about the initial momentum po distribution is lost, and the positions 
(or field amplitudes) at different times commute, 

y(k,n), y(k,T 2 ) 

The last result defines what is known as a quantum non-demolition (QND) variable, which means that one can perform 
successive measurements of this variable with arbitrary precision without modifying the wave function. Note that 
y = 8(f) is the amplitude of fluctuations of the Higgs field after inflation, so what we have found is: first, that the 
amplitude is distributed as a classical gaussian random field with probability (|31[) : and second that we can measure 
its amplitude at any time, and as much as we like, without modifying the distribution function. 

In a sense, this problem is similar to that of a free non-relativistic quantum particle, described initially by a minimum 
wave packet, with initial expectation values (x)q — xq and (p)q = po, which becomes broader by its unitary evolution, 

and at late times (t 3> mxo/po) this gaussian state becomes an exact WKB state, ^(x) = Q R X ^ 2 exp(— £1 x 2 /2), with 
Imi7 Re Q (i.e. high squeezing limit). In that limit, [x, p] « 0, and we have lost information about the initial 
position xq (instead of the initial momentum like in our case), x(t) — > p(t)t/m = pot/m and p(t) = pq. Therefore, 
not only [p(ti), pfo)] = 0, but also, at late times, [x(t\), 2(^2)] ~ 0. 



— e 

2 



cos 



(A17) 



APPENDIX B: THE WIGNER FUNCTION 



The Wigner function is the best candidate for a probability density of a quantum mechanical system in phase-space 
|39|| . Of course, we know from Quantum Mechanics that such a probability distribution function cannot exist, but the 
Wigner function is just a good approximation to that distribution. Furthermore, in the case of a gaussian state, this 
function is positive definite, and can in fact play the role of a classical probability distribution for the quantum state. 

Consider a quantum state described by a density matrix p. Then the Wigner function can be written as 



W(ylyT,plpV)= I I d fl^l e -^^^ (y-^r p y+^r) . (Bl) 



(2tt) 2 V 2 



2 



If we substitute for the state our vacuum initial condition p = \tyo)( i &o\, with \&o given by the gaussian wave function 
H30J1 . we can perform the integration explicitly to obtain 



1 ( \y 



2 



W Q {y u k ,yT ,pt,pt~) = ^2 eX H~& _4 ' /fcl ^~ ) =^>P^^'P^ 

$( yi , P = -exp(-(^-M|/ fc | 2 pf)l - Pi=Pi-T%Vi- (B2) 
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However, at time r = r , we have yl = ^7= = |/fe(r )| , p\ = J \ = (2|/ fe (r )|) \ and F fe (r ) = 0, so that p\ = p\, 
and therefore Wo describes a symmetric gaussian in phase space, with the same dispersion in both y and p-directions. 
The 2er contours of this distribution satisfy 



|/ ^+4|/ fc | 2 p1<l ^ f + l^ 1 ' i0VT = T °> (B3) 



which is a circle in phase space. On the other hand, for time t > rg, we have 



l/fcl -> — = e rk ~ yl e rk , growing mode, (B4) 

V 2k 



1 k 

^TTT ~* V o e ~ Tk ~ Pk e ~ r " > decaying mode , (B5) 

2|/fc| V 2 

so that the ellipse (|B3ll becomes highly "squeezed" , see Fig. El Note that Liouville's theorem implies that the 
volume of phase space is conserved under Hamiltonian (unitary) evolution, so that the area within the ellipse should 
be conserved, and as a consequence there is no entropy production in this process. As the probability distribution 
compresses (squeezes) along the p-direction, it expands along the y-direction. At late times, the Wigner function is 
highly concentrated around the region 

M»-isp»)'<4isr (B6) 
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FIG. 14: The 2a contour of the Wigner function l|B2fl for the mode k = 1, at times r = 0,1,2,3. It is clear that, as 
time progresses, the ellipse (IB3II becomes more elongated (squeezed), without changing its area, while the main axis rotates 
counter-clockwise. 



We can thus take the above squeezing limit in the Wigner function l|B2|) and write the exponential term as a Dirac 
delta function, 



W (y, P ) r ^°° 



1 



exp ■ 
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In this limit we have 
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(B8) 



so we recover the previous result (|A16() . This explains why we can treat the system as a classical gaussian random 
field: the amplitude of the field y is uncertain with probability distribution \dlti . but once a measurement of y is 
performed, we can automatically assign to it a definite value of the momentum, according to (|A16|) . 

Note that the condition F% 3> 1 is actually a condition between operators and their commutators/anticommutators. 
The Heisenberg uncertainty principle states that 



A^AA^B > 



!*P, B]\9) 



for any two hermitian operators (observables) in the Hilbert space of the wave function In our case, and in Fourier 
space, this corresponds to 



A^ 2 (fc)A* P 2 (fc) = F fe 2 (r) + i>i 



with |\&) = |0, tq) the vacuum wave function. On the other hand, Fk can be written as 



F k = ~ <*|p(k, t) yt (k, t) + y(k, r) pt ( k; T ) |v£r) = _ * {gk f* -f k g* k )=I m (f* g k ) , 



(B9) 
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where we have used (|20|) and a(k, ro)|^) = 0, Vk. The above relation just indicates that, for any state 'J, the 
condition of classicality (F k 3> 1) is satisfied whenever, for that state, 

(iHr), p{(r)})^{\Mr), $(r)]|)=ft=l. (Bll) 

It is this condition which allows one to substitute quantum averages of arbitrary functions G of the position and 
momentum operators by classical ensemble averages of the same function G, weighted with the Wigner probability 
distribution function, or schematically, 

<*|G(y fe , p k )\V) = J dy k dp k G{y k ,p k )Wo(y k ,Pk) + 0(h) (B12) 

yj_ 

-LJdy k G(y k ,Wy k )e I AMP , (B13) 

where we have used Eq. i|B7|l . As long as F k (r) ^> 1, we can describe the evolution of our quantum system as that 
of a classical gaussian random field. Note that, in this limit, we can ignore the normal ordering of the operators in 
G(y k , p k ). 



APPENDIX C: NON-LINEAR EVOLUTION IN PERTURBATION THEORY 

In this appendix we will give details of how to perform perturbative calculations of the non-linear evolution of 
our quantum system. We will also illustrate the perturbative expansion of the correlation functions of a classical 
random field. To simplify the expressions we will consider the case of a single component real scalar(Higgs) field, but 
generalisation to complex or multiple component case is straightforward. 

Our goal is to compute the expectation values of products of fields at different points: 

Mn,xi)...#7-„,x„)) (Ci) 

(In this section we will use the symbol 4> instead of y for the Higgs field) . Here <P(t, x) denotes the Heisenberg picture 
field operator, whose relation to the Schrodinger picture one </> s (x) is as follows: 

0(r,x)=Wt( T )^(x)W(r) (C2) 

where U (r) is the evolution operator, satisfying: 

U'{t) = -iRU{t) (C3) 

where the prime stands for derivative with respect to r and H is the full Hamiltonian. Notice that since the Hamil- 
tonian depends explicitly on time, the evolution operator cannot be written as exp{— irH}. If we set A to zero we 
get the quadratic Hamiltonian Ho considered in the gaussian approximation. The corresponding evolution operator 
is Uq(t). Now we go over to the interaction representation by writing: 

U(r)=U {r)n(T) (C4) 

where H(r) is the characteristic Moller type operator which satisfies the equation: 

n'(r) - -il4(r) H int W (r) fl = -<h£>(t) ft (C5) 

where h|J||.(t) is the interaction hamiltonian in the interaction representation. The equation for ft can be solved in 
terms of the time-ordered exponential: 

O(r) = Texpl-z fdtH^)} (C6) 
Jo 

This can be used to express the Heisenberg representation fields in terms of the (gaussian) interaction representation 
fields: 

^(r,x)=T'exp{i/ dtU^ (t)} o (r,x) Texp{-i / dtU^ (t)} (C7) 
Jo Jo 
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In the T-exponential time grows from right to left and in the T'-exponential left to right. To obtain the perturbative 
expansion one has to expand the T-exponential and H-°j. in powers of A. The latter has the form 

<iW = (^o(r,x) - 3( ^ 2T) fex)) + 0(A 2 ) (C8) 

where the second piece is the counter-term needed to renormalisc to this order. Then, substituting the expression of the 
field inside the expectation values, everything reduces to expectation values of products of interaction representation 
fields 4>o(t,x). The latter reduce, by Wick's theorem, to products of two-point functions: 

/ r ]3 h. 
e *(«-«) /*(,-)/*(/). (C9) 

This can be decomposed into a real and imaginary part. The real part corresponds to the expectation value of the 
symmetrised product, which in the gaussian theory was chosen to match with the correlation function of the classical 
random field. The imaginary part is proportional to the commutator of the fields, which is a c-number. 
We can illustrate the procedure by computing the 2-point function 

/d 3 k 
__ e ik(x-x) G(fc)Ty) (C10) 

to order A. Substituting the expression of the Heisenberg field for the t > t' case we get 

(T'exp{i£dtnW(t)} o (x, r) Texp{-i£di (t)} <fo(y> , r')T exp{-i £dtH^(t)}) 
In case r' > r the factor sitting between the two fields has to be replaced by: 

r' 

T'exp{*y dtH<J>(t)} (Cll) 

Notice the peculiar time-ordering of the operators which differs from the customary perturbative evaluation of (Fcyn- 
man) Green functions, which are time-ordered products of field operators. 

To do the calculation to order A it is better to start by expressing the Heisenberg field to this order: 

0(t,x) - <Mr,x) +J J d S [H^( S ),0o(r,x)] + ... 

= 0o (t, x) + iX J ds d 3 z [Ms, z), 0o(r, x)] (cpl(s, z) - 3(5i + S 2 s)<j)o(s, z)) (C12) 

Finally one obtains 

G(k,r,r') = f k (r)r k (r') + 6X f ds A ren (s) Im [/ fc (r)/ fc »] / fc («)/ fc *(/) 

Jo 

+6A [ T dsA tea (s) Im [f k (r')fl( S )} f k (r)fi( 8 ) . (C13) 
Jo 

The meaning of A ren (r) is given in the main text, where one can also find the symmetrised Green function. Repeating 
the calculation for a Higgs field with N c real components one gets the same expression for each component replacing 
6 by 2(iV c + 2). 

We can compare with the classical evolution. We will use the same symbol for the classical field </>(t, x). The 
equations of motion 

0" (r, x) - A0(r, x) + t0(t, x) - A0 3 (r, x) (C14) 

can be solved in perturbation theory in A. The expansion is given in terms of tree graphs with lines associated to the 
retarded propagator: 

/d 3 k 
— ^-x)l m {f k (T).n(T>)}. (C15) 
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If we now take 4>(t = 0,x), and tt(t = 0,x), to be gaussian random fields, then the field at any other time becomes 
a non-gaussian random field. The correlation functions of this field can be computed in perturbation theory by 
combining the aforementioned expansion involving the retarded propagator and the expectation value of gaussian 
random field. To match with the quantum calculation at zero-order in A, this has to be taken as the symmetrised 
version of Eq. IC9I (this is just given by the Fourier transform of the symmetric part of £): 



Notice that terms in the expansion can be associated with Feynman-type graphs, with modified rules involving two 
propagators {G^J UBS and G rct ). These, up to factors (including Heaviside 9) coincide with the real and imaginary 
parts of the quantum propagator. Indeed, the calculation of the two-point correlation function to order A matches 
exactly with the symmetrised quantum two-point function to this order. For that one has to apply exactly the same 
renormalization to the classical and quantum theories. 

Differences can arise to higher order. Essentially, the retarded- imaginary propagator in the classical theory cannot 
form loops by itself, since it arose from the expansion of the field equations. This need not be the case in the quantum 
theory. For example to second order in A there is a contribution to the two-point function given by the sunset diagram, 
with three imaginary propagators joining the two vertices. However, for low-momenta flowing through the lines (and 
large enough times) the dramatic difference in size of the real and imaginary parts of fk( T )fk( T ') justifies that the 
classical approximation would still be reasonably good. A more thorough investigation of these matters is interesting 
but exceeds the realm of this paper. 
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